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ABSTRACT 



We announce the discovery of a new low-mass, pre-main sequence eclipsing binary, MML 53. Previous observations of MML 53 
found it to be a pre-main sequence spectroscopic multiple associated with the 15-22 Myr Upper Centaurus Lupus cluster. We identify 
the object as an eclipsing binary for the first time through the analysis of multiple seasons of time series photometry from the 
SuperWASP transiting planet survey. Re-analysis of a single archive spectrum shows MML 53 to be a spatially unresolved triple 
system of young stars which all exhibit significant lithium absorption. Two of the components comprise an eclipsing binary with 
period, P = 2.097891(6) ± 0.000005 and mass ratio, q ~ 0.8. Here, we present the analysis of the discovery data. 
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1. Introduction 



stars: pre-main sequence - stars: binaries: spectroscopic - stars: fundamental parameters 



Eclipsing binaries are systems in which the masses, radii, 
temperatures, and absolute luminosities can be measured for 
two s tars with the sam e age and metallicity (Ander sen. J.L 
ll99ll:lYoung et all l200lh . In addition, apsidal motion measure- 
ments, if available, can f urther constrain interior st ellar physics 
dYoung & Arnettl l2Q05t ISchwarzschild. M.L I1957I) . Therefore, 
these objects are effective tools used extensively for calibrating 
• • i and testing theoretical stellar evolution models. 
. ' Pre-main sequence eclipsing binaries are particularly impor- 
. tant. Young stars are changing rapidly as they contract onto the 
j_i ■ main sequence and heat up, so an individual solar-type star cov- 
C$ ' ers a large range of temperatures and radii in its first 50 Myr. 
Furthermore, young stars often have star spots, coronal activity, 
and/or circumstellar dust which complicates the physics of their 
formation and early evolution. Thus, empirical measurements of 
fundamental properties of young stars with a range of masses 
and ages are needed to fully constrain the pre-main sequence 
evolution and understand how the properties of young stars are 
affected by other parameters (e.g. magnetic fields, dust, metal- 
licity, planets, multiplicity). 

Pre-main sequence eclipsing binaries are the only ob- 
jects which can provide precise measurements of the 
most fundamental physical properties of young stars, how- 
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ever these objects are extremely rare. There are only 6 
known low-mass pre-main sequence ecl ipsing binaries with 
M < 1.5 M : RXJ 0529.4+0041A dCovino et all I2000L 
12001 . VI 174 Ori (IStassun et all l 200l 2MASS J05352184- 
20130 546085 (IStassun et all 120061120071). IW 380 dlrwin et all 
120071) . Par 1802 (Cargil e et all 120081: IStassun et all 120081) . and 
ASAS J 05282 1+033 8 5 dStempels et all 120081 I n addition. 
EK Cep (Pop per. D.M.lll987|). TY CrA dCasev et afl ll998). and 
RS Cha dAlecian et al.. 2005, 2007) are known to be higher mass 
eclipsing systems with at least one PMS component. 

In this paper, we announce the discovery of a new low-mass 
PMS eclipsing binary, and the first such object discovered out- 
side of the Orion star forming region. We describe the observa- 
tions that were used to identify the object as a PMS eclipsing bi- 
nary (§0. Using the discovery data, we determine initial orbital 
parameters for the system and approximate physical properties 
of its component stars Finally, we discuss future observa- 
tions needed to fully analyse this important system (§|4). 

2. Observations 

TYC 7310-503-1 (a =14:58:37.7, 6 =-35:40:30.4), MML 53, is 
a late-type star spatially located in the region of the Scorpius 
Centarus OB association complex. It was identifed as an X- 
ray source by ROSAT (RXJ1458.6-3541) which lead to further 
spectroscopic study by several groups with the aim of identi- 
fying young stars and mapping the overall star forming region 
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Table 1 Properties of MML 53 obtained from the literature. 



Parameter 


Value 


Reference 


RA(J2000) 


14-58-37 70 




Dec(J2000) 


-35:40:30.4 




VT 


10.88 mag 


1 


I 


9 85 + 04 map 


2 


j 


O.UJ7 IT_ VJ.VJZ.'-r llldg 


3 


H 


c 062 + 055 map 


3 


K 


7.870 ± 0.026 mag 


3 


flRA 


-22.3 ±1.1 masyr 1 


4 


^DEC 


-25.0 ± 1.1 masyr 1 


4 


SpT 


K2IVe 


5 


71 


7.36 ± 0.77 mas 


5 


7 


2.0 ±3.1 km s- [ 


6 


v sin f 


30.8 ±3.1 km s" 1 


6 



(1) 

(2) 
(3) 
(4) 
(5) 
(6) 



TYCH O dH0g et all l2000h 

DENIS jThe Denis Consortiumll2005T) 
2MASS (S krutskieetalll2006h 
NOMAD jZacharias et al.Ll2004h 
dMamaiek et alll2002h~ 
SACY dTorres et all 12006) 



dWichmann et all 11997111 iMamaiek et all 120021: ITorres et all 
2006- IWhite et all l2007h . The object was first identified as a 
T Tauri star (K3-type) through the d etection of Li I /16708 ab- 
sorption by IWichmann et all d!997bl) . Ha emission and signif- 
icant lithium absorption were detected in all subsequent spec- 
troscopic observations, therefore MML 53 is confirmed to be a 
young, pre-main sequence (PMS) object. 

IMamaiek et alJ (120021) also kinematically and spatially de- 
fined MML 53 as a member of the 15-22 Myr old Upper 
Centaurus Lupus (UCL) sub-association with a 93% probabil- 
ity. The authors classified the object as a K2IVe star based on 
several spectroscopic indices and determined a kinematic paral- 
lax of 7.36 + 0.77 mas (136^3 pc). These values and the high 
probability of membership in UCL are confirmed in more re- 
cent work b y these authors (E. Mamajek priv communication). 
ITorres et"al] (2006) made a rotation measurement of vsinj = 
30.8 ±3.1 km s'and were the first to identify the target as a pos- 
sible spectrosco pic triple system w ith a systemic radial velocity, 
y — 2.0 km s" 1 . IWhite et al.l d2007l) detected a double-lined na- 
ture for the object with a single observation and measured the 
equivalent width of lithium in the individual primary and sec- 
ondary components (122 and 222 mA). I White et al.l d2007l) do 
not explicitely state which binary component is associated with 
the reported Li I EW values since they are unable to derive an 
orbit from their spectrum. The approximate values we derive for 
Li EW in Sect. l3.3l are consistent with these measurements if we 
assume the primary star has the larger EW (222 mA). 

As a known PMS object, MML 53 has also been included in 
the Spitzer Legacy project investig ating the formation and evo- 
lution of planetary systems (FEPS) iMever et all d2006l) . Spitzer 
observations of the obje ct were obtained in all IRAC and MIPS 
bands from 3.6-70 jum dCarpenter et all | 2008|). but no primor- 
dial or a debris disk was detected (Silversto ne et all 12006; 
ICarpenter et alll2009l) . In Table Q] we list the known properties 
of the star obtained from the literature. 
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Fig. 1 WASP-South photometry of MML 53 from 2006-2008 ob- 
tained in a broad V+R-band filter. The data are phase-folded 
with the ephemeris derived from the eclipse model, HJD = 
2454301.3252(7) + 0.002 + 2.097891(6) + 0.000005. 



2.1. SuperWASP Photometry 

In addition to the targeted surveys of young stars mentioned 
above, MML 53 was observed in th e field-of-view of th e 
SuperWASP transiting planet survey dPollacco et all [2006). 
SuperWASP is a wide-field photometric variability survey de- 
signed to detect transiting gas-giant planets around bright main 
sequence stars. The survey cameras repeatedly (~ 8 minute sam- 
pling) observe bright stars (V ~ 9 - 13) at high precision (1%) 
using a broad V+R-band filter. 

MML 53 wa s observed with the WASP-South telesco pe and 
instrumentation dPollacco et all 120061: IWilson et all 12008) in the 
spring observing season from 2006 to 2008. 3492 photomet- 
ric data points were obtained between 4 May - 30 July 2006; 
5602 measurements were made between 16 Feb 2007 - 19 July 
2007; and 2886 observations were taken from 18 Feb 2008 - 17 
Apr 2008. All data sets were processed independently with the 
standard SuperWASP data re duction and photometry pipeline 
(ICollier Cameron et all 12006). The light curves, which were de- 
trende d along with other s tars in the field using the SysRem algo- 
rithm dTamuz et alll2005l) . were then run through our implemen- 
tation of the box least squares (BL S) algorithm dKovacs et all 
I2002L ICoIUct Camer on et al ., 2006). The BLS algorithm is de- 
signed specifically to detect square shaped dips in brightness in 
an otherwise flat light curv e. It is very effe ctive at detecting tran- 
sits by e xtrasolar plane ts ( Tinglevl 1200 3l) and eclipsing binary 



systems (lHartman et al ., 2009)^ 

Using the BLS algorithm, we identified MML 53 as hav- 
ing an eclipsing binary light curve with a period of ~ 2.1 days. 
FigureQ]shows the phase-folded SuperWASP data. The target is 
isolated, and the field-of-view shows no other nearby stars that 
could have contaminated the large photometry aperture (~ 48") 
causing the resulting light curve. 

The primary and secondary eclipses are apparent in the light 
curve as well as sinusoidal out-of-eclipse variability which is 
particularly evident in the 2006 data. Brightness variations due 
to photospheric star spots are often present on young and/or 
active stars and can be used to measure stellar rotation peri- 
ods. Thus, we investigated the presence of rotational variabil- 
ity on MML 53 to determine if the system is tidally synchro- 
nized. We measured the period of the out-of-eclipse variability 
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Fig. 2 Lomb-Scargle periodogram results from the 2006 
SuperWASP photometry with the in-eclipse data removed. 
Periods between 0.2-30 days were searched. A periodic signal 
was detected at P ~ 2.09 days in all 3 seasons of data, but 
the signal is strongest in the 2006 data shown here. The dotted 
line denotes the orbital period of the binary, P = 2.097891(6) + 
0.000005. 



by applying the Lomb-S cargle algorithm (Sca rgle et al.L ll982; 
iHorne & Baliunasl Il986l) to each season's light curve indepen- 
dently after removing all the in-eclipse data. We searched peri- 
ods between 0.2-30 days. In the 2006 data, there is a highly sign- 
ficant sinusoidal signal with a period, P = 2.09 days (see Fig. |2j. 
Aliased peaks due to the window function are also present at a 
reduced amplitude. In the 2007 and 2008 data sets, the signal is 
weaker, but still present. This period is matched to the orbital 
period of the binary suggesting the out-of-eclipse variability is 
due to starspots on the surface of one or both of the binary com- 
ponents and that the component stars are rotating synchronously 
with the orbit. The sinusoidal variability is changing phase and 
amplitude on the timescale of months which manifests as visible 
scatter in the light curve (see Fig.Q} indicating a slow drift in the 
starspot pattern, but the measured rotation period is consistent in 
all the seasons of data. For a short period binary like MML 53, 
we expect the c omponents to be rotating synchronous ly despite 
their young age (IZahn & Bouchell ll989: MathieuL 119941 and ref- 
erences therein). 

2.2. Archive FEROS Spectrum 

After determining the system was an eclipsing binary, we 
searched for existing archival spectroscopic data on the ob- 
ject and found one high resolution (R ~ 50 000) spectrum 
located in the European Southern Observatory (ESO) archive. 
The spectrum was obtained at heliocentric Julian date, HJD = 
2453909.622300 with the FEROS spectrograph on the 2.2m 
ESO telescope at La Silla (observing program Q77 .C-0138(A), 
which forms the basis of the SACY survey, see iTorres et"al] 
(2006)). We reduced this spectrum with the echelle data reduc- 
tion package REDUCE (Pis kunov & Valeria , 120021) . using cal- 
ibration data obtained on the same night. We also downloaded 
and reduced the spectrum of a radial velocity standard, HD 
10700. Our final reduced spectrum covers the wavelength range 
of 3765-8862A and has a signal-to-noise of ~ 30 per pixel at 
6000A. 




Radial velocity (km s ) 



Fig. 3 Cross correlation function of the FEROS spectrum of 
MML 53 obtained on 2006 June 23 compared to the template. 
The analysis shows three different spatially unresolved compo- 
nents in the system with radial velocities of -85.8, 111.1 and 
-3.5 km s'for the primary, secondary, and tertiary components, 
respectively. 



When the spectrum was examined by eye, multiple stel- 
lar components were clearly present as previously reported 
(ITorres et all 120061: 1 White et all 120071) and as expected for an 
eclipsing binary wher e both the pr i mary and secondary eclipses 
are visible. However, I White etaD (120071) report t he system as 
a dou ble-lined spectroscopic binary (SB2) while ITorres etaD, 
(2006) suggest the possibility of a third component (SB3). To 
determine whether the system contains a third star and to mea- 
sure the radial velocities of the components, we perform a stan- 
dard cross-correlation (Fig. [3} and also apply our implementa- 
tion of the least-squares deconvolution (LSD) algorithm to the 
MML 53 spectrum. Least-squares deconvolution allows for de- 
riving a very high signal-to-noise average absorption line pro- 
file of the system by properly combi ning the many indi vidual 
lines throughout the entire spectrum (Donati et al., 1997). This 
technique can be used in spectroscopic binary work to identify 
faint unresolved companions. Both the cross-correlation func- 
tion (CCF) and the result of the LSD analysis clearly show three 
stellar components present in the spectrum. We measure their 
radial velocities to be -85.8, 111.1 and -3.5 km s _1 , for the pri- 
mary, secondary and tertiary, respectively. 



3. Analysis 

3. 1 . Light Curve Modelling 

In order to derive an ephemeris for the system and initial parame- 
ters for the individual eclipsing binary components, we first gen- 
erated a single recitified light curve from the discovery data by 
removing the out-of-eclipse variability on each night and com- 
bining all 3 seasons of SuperWASP photometry. We fit a sec- 
ond order polynomial to the out-of-eclipse photometry on each 
night and applied the polynomial fit to all the data obtained that 
night. We excluded from the rectified light curve any data taken 
on nights in which there was no out-of-eclipse photometry. The 
fitting functions are not physical, and we made no attempt to 
model the starspot variability because the existing data are only 
in a single band and therefore any physical model would be too 
full of degeneracies to provide a useful result. The data include 
17 near complete primary or secondary eclipses which show the 
eclipse minimum and out-of-eclipse photometry before and/or 
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The model fit also gives a measurement for the relative 




sum of the radii, (Ri + R2)/a = 0.260, the inclination angle, 
i - 83.1°, and the surface brightness ratio in the SuperWASP 
filter, Jv+r = 0.461. These values depend on the flux contribu- 
tion of the third component which we estimate to be 15% of the 
total luminosity based on our self consistent analysis of the light 
curve and the FEROS spectrum described in Sect. 13.31 However, 
we stress that the fitted parameters will change when definitive 
temperatures for all three components and accurate luminosity 
ratios can be derived through spectral disentangling of multi- 
ple high resolution, high signal-to-noise spectra of MML 53. 
Furthermore, star spots which we know are present are likely 
to have an affect on these parameters, and multi-band photome- 
try is necessary to derive a comprehensive solution that models 
both the spots and the eclipses. 



phase 



Fig. 4 Rectified SuperWASP light curve phase-folded with 
ephemeris, HJD = 2454301.3252(7) + 0.002 + 2.097891(6) + 
0.000005. The EBOP model fit is overplotted (solid black line). 

after the eclipse. The final rectified light curve contains 11808 
photometric measurements. 

We fit the rectified light curve u sing the JKT E cl ipsing 
Binary Or bit Program (EBOP) ( Pop per & Etzel 1 198 it 
ISouthworthl |2007l) . The program determines the optimal 
model light curve that matches the observed photometry and 
reports the binary parameters for the model. The algorithm on 
which it is based is only valid when analysing well-detached 
eclipsing binaries in which the tidal distortion is small (i.e. 
nearly spherical stars with oblateness < 0.04). This is the 
case for MML 53. The derived light curve parameters include 
the period, P, time of minimum light, To, surface brightness 
ratio in the SuperWASP filter, Jv+r, relative sum of the radii, 
(Ri + Rz)/a, inclination angle, i, eccentricity, e, and angle of 
periastron, u. The routine takes into account the effects of limb 
darkening, gravity brightening, and reflection effects and can 
account for the presence of light from a third component. 

We adopt ed the quadratic limb darkening coefficients from 
IClared <l2000h using the temperatures for the eclipsing compo- 
nents that are defined below in Sect. 13.31 r e ff,i = 4886 K and 
T e s,2 = 4309 K. We initially ran the fitting program allowing the 
eccentricity to be a free parameter, however the resulting value 
was within 3cr of zero, and there is no further evidence in the 
light curve for a non-circular orbit. 

Thus, in the final run of the light curve modelling program, 
we fixed the eccentricity and angle of periastron to zero. The 
result is shown in Fig. [4] with the best fitting model overplotted 
on the final rectified light curve. The fit provides a precise system 
ephemeris of: 

Min(HJD) = 2454301. 3252(7)+0.002+2.097891(6)±0.000005£ 

As an additional check, we determined the ephemeris from 
each individual season of data, and found the periods agree to 
within 5 x 10~ 6 days. However, the epoch of mininum light varies 
by ~ 0.001 in phase (~ 3 min) between the 2006 and 2007-08 
seasons. This offset could be caused by effects from the third 
component (e.g. direct gravitational influence, light travel time 
variations), if ultimately confirmed with current epoch eclipse 
data. Based on the variations from season to season, we adopt an 
uncertainty on the time of minimum light of 0.002 days and an 
uncertainty on the period of 0.000005 days. 



3.2. Initial Estimates of the Masses and Radii of the MML 53 
Eclipsing Components 

Using the radial velocity measurements of the primary and 
secondary star from the single archive spectrum, the precise 
ephemeris and inclination angle estimate from the photome- 
try, and the systemtic radial velocity from the literature, we 
derive approximate masses for the individual eclipsing com- 
ponents of MML 53. The FEROS observation was obtained 
very near to quadrature at an orbital phase of 0.287. With the 
phase of the observation fixed, the radial velocity measure- 
ments for the primary and secondary (-85.8, +111.1 km s _1 ) 
constrain the amplitudes of the sinusoidal (circular orbit) radial 
velocity cu rves for the two c omponents, K\ and K 2 . Adopting 
2.0 km s' dTorres et all 120061) for the systemic RV, we find 
K\ = 90 and K2 = 112 km s _1 (assuming no uncertainty in 
the RV measurements). This gives a mass ratio for the system, 
M%jM\ = 0.8. Taking the period and inclination angle derived 
from the SuperWASP eclipse photometry, we find individual 
masses of M\ ~ 1.0 M , M 2 ~ 0.8 M , and an orbital sepa- 
ration, a ~ 8.4 R Q . 

The derived mass estimates are consistent with the pre-main 
sequence status of MML 53 determined from the observed youth 
indicators previously discussed. Young stars should be compar- 
atively large in radius, since they are still contracting onto the 
main sequence. For MML 53, the sum of the radii measured 
from the light curve (Ri + R2 = 2.2 Rq) is ~ 30% larger than 
what is expected for two main sequence stars (age of 300 Myr) 
with masses of 1. and 0.8 Mq accor ding to theoretical stellar 
evolution models (Baiaff eetailll998h . 

Finally, we note the properties of MML 53 derived here are 
only approximations, and a complete radial velocity curve will 
be required to determine accurate masses for the stellar com- 
ponents of MML 53. Furthermore, the third component, which 
cannot be studied in detail with the existing data could also affect 
the radial velocity measurements if it is found to be gravitation- 
ally bound to the system, which is still unknown. 

3.3. Composite Model Spectrum: Temperatures and Lithium 
Abundances 

Despite the lack of time-series spectra, we perform a joint analy- 
sis of the single observed FEROS spectrum and the SuperWASP 
light curve in order to approximate the effective temperatures, 
radii and relative luminosities of the three stellar components of 
MML 53. To do this, w e derive a three -component model spec- 
trum (based on synthetic Kurucz ( 1993) model atmospheres) that 
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Wavelength (A) 

Fig. 5 FEROS spectrum of the order containing Ha and the Li i 
doublet at 6708A. The solid black line is the observed spectrum. 
Overplotted is the three-component model spectrum (grey line). 
The model is a good match to the observed spectrum except near 
Ha where MML 53 shows emission and Li i where it shows sig- 
nificant absorption. Plotted below are the three individual syn- 
thetic spectra with temperatures of 4886 K, 4309 K, and 4130 K 
for the primary, secondary and tertiary. They are scaled accord- 
ing to the derived R-band luminosity ratios of 0.64:0.21 :0. 15 and 
shifted by the velocity offsets found in the CCF analysis. 



is a reasonable fit to the observed spectrum and is also consis- 
tent with the constraints from the light curve. The model for each 
component is defined by its effective temperature and its lumi- 
nosity relative to the total system luminosity at ~ 6700A (R- 
band). These properties are interdependent, so we iterate until 
they provide a self consistent solution. 

After several iterations, we find temperatures of 4886 K, 
4309 K, and 4130 K for the primary, secondary and tertiary, 
respectively. In addition, the fractional R-band luminosities we 
use to scale the model spectra are 0.64:0.21:0.15 for the three 
components. These values provide a self consistent solution in 
which the third light input value for the light curve model is 
reproduced by the final temperature and radius of the tertiary. 
In addition, the properties of the eclipsing components are con- 
sistent with the light curve model. Finally, these values result 
in a composite model spectrum which is a good match to the 
absorption line depths of all three components in the observed 
spectrum in the region between Ha and Li i. In Fig. [5] we show 
this region of the observed spectrum (black solid line) overplot- 
ted with the composite model defined by our best fitting stellar 
properties. The model assumes a solar metallicity for all three 
components, and we adopt rotational velocities for the stars of 
vsin/ = 29.0, 23.2, and 20.0 km s _1 . Note that the vsin/ val- 
ues for the primary and the secondary are determined from the 
assumption of corotation and are reasonably consiste nt with the 
value of 30.8 km s 'measured bv lTorres et al.l (|2006). 

As can be seen in Fig. [5] MML 53 shows significant Li i 
absorption in all three components and Ha in emission for the 
secondary and tertiary. We measure the equivalent width (EW) 
for each Li line in the combined spectrum to be 236, 75, and 
82+ 10 mA for the primary, secondary and tertiary, respectively. 
We then correct each of these values for the increased con- 
tiuum due to the other components to derive corrected equivalent 



widths of EW(Li)= 369+ 15, 356 + 47, 550 ± 67 mA which we 
translate into lithium abundances of log n (Li) = 3.2, 2.3, and 
3.0. 

The uncertainties on the lithium EWs quoted above do not 
include the systematic error in the relative luminosity ratios 
which are difficult to estimate. The lithium abundance of the 
tertiary is unexpectedly larger than that of the secondary which 
suggests some inaccuracies in our derived temperatures and lu- 
minosities. However, with additional observations of time series 
spectra, we will be able to apply spectral disentangling which 
will allow for deriving accurate temperatures for the individual 
components and measuring Li equivalent widths to greater pre- 
cision. 



4. Summary 

We report on the discovery of a new pre-main sequence ecips- 
ing binary, MML 53, located in the Lupus star forming region. 
There are very few such objects known, but they are extremely 
important for informing and calibrating stellar evolution models 
at young ages where stars are changing rapidly as they contract 
onto the main sequence. This is only the seventh low-mass PMS 
eclipsing binary to be discovered. 

Previous studies found the object to be a young spec- 
troscopic multiple system associated with the 15-22 Myr 
Upper Centaurus Lupus cluster. Through our analysis of 
SuperWASP photomery and a single archive spectrum, we de- 
termine MML 53 to be a triple system of young stars where two 
of the components are eclipsing. All three components show sig- 
nificant lithium absorption, and Ha emission is seen in the sec- 
ondary and tertiary. 

We determine the orbital period of the eclipsing pair to be 
P = 2.097891(6) + 0.000005 using many primary and secondary 
eclipses observed with the SuperWASP instrument over several 
years. The light curve also shows out-of-eclipse variability con- 
sistent with star spots on the surface of at least one of the eclips- 
ing components which is rotating synchronously with the orbit. 

Analysis of the existing data suggests the eclipsing compo- 
nent stars have masses of Mi ~ 1.0 M Q and M 2 ~ 0.8 M Q . 
Therefore, MML 53 is a slightly older analogue of the VI 174 Ori 
PMS eclipsing binary found in the ~ 10 Myr Orion OB lc sub as- 
sociation (Mi = 1.0 M and M 2 = 0.73 M Q ). When MML 53 is 
fully analysed, relative comparison of the two objects will allow 
for constraining the PMS evolution of a solar-type star to very 
high accuracy. In addition, because MML 53 has an expected 
age older than any of the other known low-mass PMS eclipsing 
binaries, it will be valuable for exploring a slightly later stage of 
PMS evolution. 

However, additional data is required to derive precise fun- 
damental properties for the components of MML 53 before 
comparing to stellar evolution models and other PMS eclips- 
ing binaries. Time-series spectroscopic data are needed to de- 
fine a complete radial velocity curve which will allow for de- 
termining masses of the primary and secondary components. 
Additional high resolution, high signal-to-noise spectra will al- 
low for determining precise temperatures and relative luminosi- 
ties of the three components stars. Multi-band photometry will 
ultimately provide a complete solution for the system with in- 
dividual masses, radii, temperatures and luminosities for two 
young, pre-main sequence stars. 
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